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LONG-TERM  GOALS 

The  long-tenn  goals  of  this  project  include  the  development  of  objective  and  efficient  tools  for 
designing  optimal  sampling  strategies  for  oceanic  applications;  and  to  develop  and  apply  practical 
methods  for  ocean  data  assimilation. 

OBJECTIVES 

The  immediate  scientific  objective  of  this  research  is  to  develop  an  ensemble -based  tool  for  the 
objective  design  of  an  optimal  sampling  strategy  for  long-term  monitoring  programs.  Subsequently, 
this  project  aims  to  apply  this  tool  for  adaptive  sampling  to  reduce  errors  in  data  assimilating 
predictive  models. 

APPROACH 

We  have  developed  a  method  for  optimal  array  design  and  applied  it  to  a  suite  of  applications, 
including  the  design  of  a  surface  mooring  array  in  the  tropical  Indian  Ocean  (Sakov  and  Oke  2007). 
The  method  builds  on  the  work  of  Bishop  et  al.  (2001),  using  data  assimilation  theory  to  determine  the 
observation  locations  that  best  constrain  a  data  assimilating  ocean  model.  The  method  seeks  to  identify 
the  set  of  observations  that  minimize  the  analysis  error  variance  of  a  predefined  variable  or  quantity. 
We  assume  that  the  purpose  of  making  a  series  of  observations  is  to  initialize  a  data  assimilating  ocean 
model.  At  present,  we  assume  that  the  data  assimilation  scheme  used  to  assimilate  the  observations  is 
ensemble  optimal  interpolation  (EnOI;  e.g.,  Oke  et  al.  2006).  Under  this  project,  we  propose  to  modify 
our  approach,  to  accommodate  other  data  assimilation  methods  (e.g.,  Ensemble  Kalman  Filter,  EnKF). 

The  method  we  describe  here  can  be  generalized  to  employ  either  an  ensemble  of  model  states  from  a 
long  model  run,  or  from  an  ensemble  of  model  runs;  or  time  series  of  observations  (e.g.,  gridded  sea- 
level  or  sea  surface  temperature). 

An  important  consideration  in  the  development  of  optimal  sampling  strategies  is  an  understanding  of 
the  impact  of  different  observation  types  on  a  data  assimilating  predictive  model.  Under  this  and  a 
related  project  (Bluelink)  we  have  explored  this  question  (see  Oke  and  Schiller  2007,  GRL). 
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Public  reporting  burden  for  the  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources,  gathering  and 
maintaining  the  data  needed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this  collection  of  information, 
including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite  1204,  Arlington 
VA  22202-4302.  Respondents  should  be  aware  that  notwithstanding  any  other  provision  of  law,  no  person  shall  be  subject  to  a  penalty  for  failing  to  comply  with  a  collection  of  information  if  it 
does  not  display  a  currently  valid  OMB  control  number. 
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WORK  COMPLETED 

We  have  applied  an  ensemble-based  method  for  optimal  array  design  to  several  applications,  including 
the  design  of  a  surface  mooring  array  in  the  tropical  Indian  Ocean  (Sakov  and  Oke  2007).  In  addition 
to  this  study,  we  have  extended  the  method  for  multivariate  and  three-dimensional  applications. 

We  have  conducted  a  series  of  observing  system  experiments  (OSEs)  to  investigate  the  relative 
importance  of  different  observation  types  (Argo,  altimeter  and  SST  data)  for  constraining  an  eddy 
resolving  ocean  model. 

RESULTS 

Objective  array  design: 

We  have  developed  a  simple,  versatile,  computationally  efficient  ensemble-based  method  for 
objectively  designing  an  observation  array.  The  versatility  of  the  method  is  demonstrated  by  Sakov 
and  Oke  (2007)  through  a  series  of  applications  to  the  tropical  Indian  Ocean  (TIO).  The  first 
application  uses  model-generated  fields  of  high-pass  filtered  mixed  layer  depth  to  design  an  array  to 
monitor  intraseasonal  variability.  The  second  uses  gridded  observations  of  sea-level  anomaly  to  design 
an  array  to  monitor  intraseasonal  to  interannual  variability.  For  both  applications,  the  objectively 
designed  arrays  are  compared  to  an  array  that  will  soon  be  implemented  under  the  auspices  of  the 
CLIVAR-GOOS  Indian  Ocean  Panel  (CG-IOP).  We  conclude  that  the  CG-IOP  array  produces  results 
that  compare  well  to  the  objectively  designed  arrays  for  intraseasonal  variability,  and  that  observations 
to  the  east  and  north-east  of  the  TIO  and  south  of  India  are  most  important  for  resolving  intraseasonal 
variability.  We  also  find  that  observations  near  9°S,  where  seasonal  Rossby  waves  dominate,  are 
important  for  observing  seasonal  to  interannual  variability.  The  described  method  for  objective  array 
design  can  be  applied  to  a  wide  range  of  geophysical  applications  where  time  series  of  gridded 
modelled  or  observed  fields  are  available. 

Observing  System  Experiments: 

We  wish  to  quantify  the  relative  importance  of  different  observation  types,  namely  Argo,  altimeter  and 
SST  data,  for  constraining  an  eddy  resolving  ocean  reanalysis.  We  use  the  Bluelink  forecast  system 
(see  related  projects,  below;  and  Oke  et  al.  2005;  2007)  to  conduct  a  series  of  OSEs  to  address  this. 
Results  from  this  study  are  presented  by  Oke  and  Schiller  (2007).  A  qualitative  assessment  of  the 
impact  of  with-holding  each  data  type  in  the  Tasman  Sea  is  shown  in  Figure  1.  The  circulation  in  the 
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Tasman  Sea  is  dominated  by  the  East  Australian  Current  (EAC)  and  a  rich  field  of  mesoscale  eddies. 

In  Figure  1,  we  show  6-day  composite  SST  from  AVHRR  measurements  (first  column)  that  are  not 
assimilated,  and  6-day  averaged  SST  for  each  OSE  (columns  2-6).  The  vector  paths  plotted  in  Figure  4 
show  Lagrangian  trajectories  over  each  6-day  period,  computed  from  the  time-dependent  reanalysed 
surface  currents.  Fields  are  shown  for  each  OSE  in  mid-January,  mid-February  and  mid-March  2006. 
The  fields  in  Figure  4  show  that  when  SST  is  assimilated,  the  reanalysed  SST  is  very  similar  to  the 
observed  SST.  By  contrast,  when  SST  is  not  assimilated  the  reanalysed  SST  differs  significantly  from 
observations.  These  differences  are  reduced  in  AFTIM+Argo,  when  altimetry  and  Argo  data  are 
assimilated,  compared  to  NONE,  when  no  data  are  assimilated.  It  is  clear  from  Figure  4  that  the 
mesoscale  circulation  is  qualitatively  similar  for  each  OSE  when  altimetry  is  assimilated,  but  is  quite 
different  otherwise.  For  example,  the  position  of  the  Tasman  front  in  January  and  February,  denoted  by 
the  separation  of  the  EAC  from  the  coast  around  32°S,  is  very  similar  in  the  OSEs  that  assimilate 
altimetry.  Also  note  the  robustness  of  the  anti-cyclonic  circulation  around  152°E  and  37°S  in  February 
and  March;  and  the  coastal  meander  near  29°S  in  February  that  is  evident  in  the  observations  and  the 
OSEs  that  assimilate  altimetry.  These  considerations  lead  us  to  the  conclusion  that  altimetry  is  required 
for  mesoscale  circulation  to  be  resolved. 
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Figure  1:  Observed  SST  from  6-day  composite  A  VHRR  (left  column);  and  5-day  averaged  SST  with 
virtual  drifter  paths  overlayed  for  each  OSE  (columns  2-6)  in  the  Tasman  Sea  for  mid-January 
(top),  mid-February  (middle)  and  mid-March  (bottom)  of 2006  (adapted  from  Oke  and  Schiller 
2007).  These  results  demonstrate  the  importance  of  altimeter  data  for  constraining  the  mesoscale 
ocean  circulation;  and  the  importance  of  SST  for  constraining  the  details  of  the  mesoscale  SST 
variability.  [Figure  showing  comparisons  between  observed  SST  and  SST  from  each  OSE] 
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IMPACT/APPLICATIONS 


The  tools  for  objective  array  designed  that  have  been  developed  under  this  project  are  potentially  very 
powerful.  Given  a  time  series,  or  model  ensemble,  of  oceanic  fields  for  some  region,  these  tools  enable 
us  to  efficiently  design  an  observation  array  for  that  region.  Our  experiments  indicate  that  the  most 
critical  aspect  of  any  application  is  formulating  the  cost  function.  That  is,  in  determining  exactly  what 
it  is  that  we  wish  to  monitor.  For  example,  an  "optimal”  array  for  monitoring  intraseasonal  mixed 
layer  depth  is  likely  to  be  very  different  from  an  "optimal”  array  for  monitoring  interannual 
variability.  This  is  demonstrated  by  Sakov  and  Oke  (2007). 

The  results  from  the  OSEs  presented  by  Oke  and  Schiller  (2007)  provide  an  important  assessment  of 
the  performance  of  the  Global  Ocean  Observing  System  (GOOS)  for  a  data-assimilating,  eddy¬ 
resolving  ocean  forecast  system.  Such  information  is  essential  for  planning  investment  in  the 
maintenance  and  development  of  the  GOOS. 

RELATED  PROJECTS 

Bluelink,  is  a  partnership  between  CSIRO,  the  Bureau  of  Meteorology  and  the  Royal  Australian  Navy. 
Many  of  the  research  activities  undertaken  in  Bluelink  have  strong  synergies  for  the  project  that  is  the 
subject  of  this  annual  report.  The  main  objective  of  Bluelink  is  the  development  and  application  of  an 
ocean  forecast  system  for  the  mesoscale  circulation  around  Australia.  This  system  is  comprised  of  a 
global  configuration  of  the  Modular  Ocean  Model  (version  4.0d;  Griffies  et  al.  2004),  with  1/10° 
resolution  around  Australia,  and  the  Bluelink  Ocean  Data  Assimilation  System  (BOD AS;  Oke  et  al. 
2005).  Briefly,  BODAS  is  an  EnOI-based  system  that  can  assimilate  observations  of  SLA  from 
satellite  altimetry  or  coastal  tide  gauge  stations,  satellite  SST,  and  in  situ  temperature  and  salinity  from 
Argo  etc.  The  details  of  BODAS  are  presented  by  Oke  et  al.  (2007),  along  with  a  comprehensive 
assessment  of  the  latest  Bluelink  ReANalysis  (BRAN)  experiment.  Oke  et  al.  (2007)  show  that  that 
reanalysed  fields  are  typically  within  6-12  cm  of  withheld  altimetric  observations,  within  0.5-0. 9°  of 
observed  sea-surface  temperature  and  within  4-7  cm  of  observed  coastal  sea-level.  Comparisons  with 
Argo  profiles  and  surface  drifting  buoys  show  that  BRAN  fields  are  within  1°  of  observed  sub-surface 
temperature,  within  0.15  psu  of  observed  sub-surface  salinity  and  within  0.2  m/s  of  near-surface 
currents.  An  example  of  the  comparisons  presented  by  Oke  et  al.  (2007)  is  shown  in  Figure  2,  where 
observed  and  reanalysed  SST  are  compared  off  Western  Australia  during  2004.  The  observations  used 
in  these  comparisons  were  not  assimilated  into  BRAN.  The  agreement  between  these  fields 
demonstrates  that  the  Bluelink  system  can  realistically  represent  the  mesoscale  ocean  circulation.  Oke 
et  al.  (2007)  present  a  more  comprehensive  assessment  of  this  system,  including  an  assessment  of  the 
short-range  forecast  skill  of  the  Bluelink  system  for  SST  and  SLA. 

The  Bluelink  Ocean  Forecasting  System  became  operational  at  the  Bureau  of  Meteorology  in  August 
2007  (www.bom.gov.au/oceanography/forecasts/).  Subsequently,  this  system  has  produced  two  7-day 
forecasts  each  week. 
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Figure  2:  Comparisons  between  observed  and  reanalyzed  SST,  demonstrating  the  good  quantitative 
skill  of  the  Bluelink  forecast  system.  Columns  1,  3  and  5  show  6-day  composite  SST  from  AVHRR 
measurements.  Column  2,  4  and  6  show  the  corresponding  5-day  averaged  SST  from  BRAN. 
Comparisons  are  shown  for  the  15th  day  of  each  month  of 2004,  as  labelled  in  the  panels  showing 
observations.  Overlaying  the  BRAN  SST  are  Lagrangian  trajectories,  computed  from  the  time- 
dependent  surface  velocities  from  BRAN  over  the  5-day  averaging  period.  Overlaying  the  observed 
SST  are  paths  of  observed  surface  drifting  buoys;  adapted  from  Oke  et  al.  (2007). 

[Figure  showing  comparisons  between  observed  and  reanalysed  SST] 
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